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Abstract. S-Shock wave type solutions in the multidimensional system of con- 
servation laws 

pt + V-(pF(U))=0, (pU)t+V-(pN(U)) = 0, 

are studied, where F = (Fj) is a given vector field, N = (Njk) is a given 
tensor field, Fj, N kj : R n -> R, j, k = 1, . . . , n; p(x,t) G K, U(x,t) G M". The 
well-known particular cases of this system are zero-pressure gas dynamics in a 
standard form 

p t + V ■ {pU) = 0, (pU) t + V • (pU ® U) = 0, 
and in the relativistic form 

Pt + V ■ { P C{U)) = 0, (pU) t + V-{pU® C{U)) = 0, 
where C(U) = , C ° U = , en is the speed of light. We introduce the integral 

identities which constitute definition of 5-shocks for the above systems and 
using this definition derive the Rankine— Hugoniot conditions for curvilinear 8- 
shocks. We show that <5-shocks are connected with transportation processes and 
concentration processes and derive the <5-shock balance laws describing mass and 
momentum transportation between the volume outside the wave front and the 
wave front. In the case of zero-pressure gas dynamics the transportation process 
is the concentration process. We also prove that energy of the volume outside 
the wave front and total energy are nonincreasing quantities. The possibility of 
the effect of kinematic self- gravitation and the effect of dimensional bifurcations 
of 8-shock in zero-pressure gas dynamics are discussed. 



1. Introduction 

1.1. L°°-type solutions. As is well known, even in the case of smooth (and, cer- 
tainly, in the case of discontinuous) initial data U°(x), we cannot in general find a 
smooth solution of the one dimensional system of conservation laws: 

f U t + (F(U)) x = 0, in Mx (0, oo), 

\ U = U°, inKx {£ = ()}, 

where F : R m -> R m is called the flux-function associated with (JTTTJ) ; U° : R -> 
M. m are given smooth vector- functions; U = U(x,t) = (ui(x,t), . . . ,u m (x,t)) is the 
unknown function, x € R, t > 0. 

Quoting from Evans' book, "the great difficulty in this subject is discovering a 
proper notion of weak solution for the initial problem (|l.ip " [T9l 11.1.1.]. "We must 
devise some way to interpret a less regular function U as somehow "solving" this 
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initial-value problem. But as it stands, the PDE does not even make sense unless U 
is diffcrentiablc. However, observe that if we temporarily assume U is smooth, we 
can as follows rewrite, so that the resulting expression does not directly involve the 
derivatives of U. The idea is to multiply the PDE in by a smooth function 

ip and then to integrate by parts, thereby transferring the derivatives onto </?" [HI 
3. 4.1. a.]. According to the above reasoning, the following definition is introduced: it 
is said that U £ L°° (ffi x (0, oo); R m ) is a generalized solution of the Cauchy problem 
(II. ip if the integral identities 

U-<p t + F(U)-!p x )dxdt + J U°{x) ■ p(x,0)dx = (1.2) 

hold for all compactly supported test vector- functions ^:Rx [0, oo) — > K m , where 
is the scalar product of vectors, and J f(x) dx denotes the improper integral 

S-co f( X ) dx - 

Using Definition (|1.2|) . one can derive the classical Rankinc- Hugoniot conditions 
for shocks (see, e.g., [HI 11.1.1.]). 

1.2. (5'"'-Shock wave type solutions, n = 0, 1, It is well known that there are 

"nonclassical" situations where, in contrast to Lax's and Glimm's classical results, 
the Cauchy problem for a system of conservation laws does not possess a weak L°° - 
solution or possesses it for some particular initial data. In order to solve the Cauchy 
problem in this "nonclassical" situation, it is necessary to introduce new singular 
solutions called S -shocks (see [1], [§]- [10], [13]- [17], [22], [24]- [31], [39]- [41], [44]- 
|49| and the references therein) , which is a solution such that its components contain 
Dirac measures. 

To illustrate the above remark, we consider the Ricmann problem for the one- 
dimensional system of zero-pressure gas dynamics 

Pt + (pu) x = 0, (pu) t + (pu 2 ) x = 0. (1.3) 

If we seek a shock solution of this problem 

u(x, t)=u + + [u]H(-x + </>(*)), p(x, t)= p+ + [p]H(-x + <t>(t)), 

we can easily verify that the velocity of the shock is 

_ W] 



where [u] = u_ — u + . [p] — p_ — p + . The last relations imply that 

P+P-{u- - u+) 2 = 0. 

Thus, in general, if u_ ^ u + the Riemann problem has no shock solutions. As was 
written in the excellent paper by A.N. Kraiko [26] page 502], to construct a solution 
for this case, it is necessary to introduce nonclassical discontinuities which may 
carry mass, momentum, energy. It is in the above-mentioned above "nonclassical" 
situations where we have to introduce i5-shocks to solve this Riemann problem for 
general initial data. Indeed, according to the above quoted papers, in the class of 
(5-shock type solutions this Ricmann problem has a solution for general initial data 
(in particular, see [HI Sec. 4]). 

The theory of <5-shocks has been intensively developed in the last ten years. 

Recently, in [36], [37], [43], a concept of <r n ) -shock wave type solutions was intro- 
duced, n = 1,2, ... . It is a new type of singular solution of a system of conservation 
laws such that its components contain delta functions and their derivatives up to 
n-th order. In [36] the theory of 5' -shocks was established. The results [36], [37], [43] 
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show that systems of conservation laws can develop not only Dirac measures (as in 
the case of (5-shocks) but their derivatives as well. 

The above singular solutions are connected with transportation processes and 
concentration processes pQ , [9] , [10] , [36] , [45] . 

8- and (5(™)-shocks, n = 1, . . . , do not satisfy the standard L°°-integral identities 
(|1.2p . Consequently, to deal with these singular solutions, we need 

• to discover a proper notion of a singular solution and to define in which sense 
it may satisfy a nonlinear system; 

• to devise some way to define a singular superposition (product) of distributions 
(for example, the product of the Hcaviside function and the delta function). 

Unfortunately, using the above cited instructions from the Evans' book [TIB 
3.4.1. a.], <5(") -shock wave type solutions cannot be defined. Indeed, if by integrating 
by parts we transfer the derivatives onto a test function ip, under the integral sign 
there still remain nonlinear terms undefined in the distributional sense, since the 
components of a solution may contain Dirac measures and their derivatives. 

Thus we need to develop a special technique. Several approaches to solving 8- 
shock problems are known (see pQ, [8]- [10], [13]- pj], [22], [24], [28]- [31], [39]- 
[41], [11]- [46], [48], [H] and the references therein). One of them was proposed 
in [Hj- [U]. In these papers the weak asymptotics method) for studying the dynam- 
ics of propagation and interaction of different singularities of quasi-linear partial 
differential equations and systems of conservation laws was developed. It appears 
that the weak asymptotics method is a proper technique to deal with 8- and (5'-shocks. 
In the framework of the weak asymptotics method, in [14j- [16j . [44] . definitions of 
(5-shock wave type solutions by integral identities for systems of conservation laws 
and in [36], [37] the corresponding definition of (5'-shock wave type solutions were 
introduced. These definitions give natural generalizations of the classical definition 
of the weak L°°-solutions (|1.2|) relevant for the structure of 8- and (5'-shocks. If a 
solution of the Cauchy problems contains no 8 and (5'-terms then these definitions 
coincide with the classical definition (|1.2I) . In |TJ, [13]- [16], [36], [37], [S], by using 
this technique, some Cauchy problems admitting 8- and ^'-shocks were solved. As 
far as we know, some problems related to 8- and ^'-shocks can be solved only by 
using the weak asymptotics method. 

In the numerous papers cited above (5-shocks in the system of zero-pressure gas 
dynamics were studied. This is related with the fact that this system has a physical 
context and is used in applications (see bellow). 

In p~7] , for the one-dimensional case of zero-pressure gas dynamics t|l ,3|) a global 
(5-shock wave type solution in the sense of Radon measures was obtained. In [22] . 
for this system the uniqueness of the weak solution is proved for the case when the 
initial value is a Radon measure. 

The multidimensional zero-pressure gas dynamics has the form 

p t + V • (pU) = 0, (pU) t +V ■( P U®U)=0, (1.4) 

where p = p(x, t) > is the density, U = (ui(x, t), . . . , u n (x, t)) € R" is the velocity, 
V = ( Tj^-, . . . , gf-) , • is the scalar product of vectors, <X> is the usual tensor product 
of vectors. 

In [29]- [31], [16], [49] the planar (5-shock wave type solution in (|1.4|) is defined as 
a measure-valued solution. The measure-valued solution is defined in the following 
way. Let BM(M. n ) be the space of bounded Borcl measures on R n . A pair (p,U), 
where p(x,t) G C(BM(R n ), [0, oo)), U(x,t) e (L°° (L°°(R n ), [0, oo)))", and U is 
measurable with respect to p at almost all t > 0, is said to be a measure-valued 
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solution of (|1.4[) in the sense of measure if 




(1.5) 



hold for all <p(x,t) G V(R n x [0, oo)). 

In this approach a smooth discontinuity surface £ is parametrized as X = X(s), 
t = t(s) (s E R n ), separating (X, i)-space into two infinite parts f2i and f^, N = 
(Nx, Nt) is the space-time normal to the surface E. The delta-shock solution takes 
the form 



(p,U)(X,t))= { (w(s,t)S(X-X(s,t),t),Ug(s,t)), (X,t) e E, (1.6) 



Here Ug is the velocity at the points of discontinuity, (pi, U\) and (p%, U-i) are 
smooth solutions of (| 1 .4|) in the regions Oi and f^2 respectively. 

In [39], [40], for the 2-D case of system (|1.4|) a notion of generalized solutions 
in terms of the Radon measures is introduced, and the problem of the propagation 
of 5-shock waves is considered. The existence of a global weak solution for the 
multidimensional system of "zero-pressure gas dynamics" is obtained in |41j . The 
approach of the latter paper is based on introducing of Lagrangian coordinates and 
on the Dafcrmos entropy condition. In pQ, for a multidimensional system of zero- 
pressure gas dynamics in non-conservative form 



the Cauchy problem related with propagation of a 5-shock wave was solved. In 1 1 1 j , 
for multidimensional continuity equation (the first equation in system (|1.4[) ) the 
possibility of existence of <5-shock was considered. 

1.3. The physical context of zero-pressure gas dynamics. Study of zero- 
pressure gas dynamics and its generalization is important for applications. The 
zero-pressure gas dynamics can be considered as a model of the "sticky particle dy- 
namics" . These models are used in many different areas of physics. Zero-pressure 
gas dynamics was used to describe the formation of large-scale structures of the 
universe [42] . |50j ; in a mathematical modeling of prcssureless mediums, in models 
of dusty gases (see the excellent papers by A.N. Kraiko and collaborators [35], [2"6]). 
in modeling two-phase flows with solid particles or droplets (see the well-known pa- 
pers by A.N. Osiptsov and collaborators [33]- [35]). The presence of particles or 
droplets may drastically modify the flow parameters. Moreover, large number of 
phenomena that are absent in pure gas flow are inherent in two-phase flows. Among 
them there are local accumulation and focusing of particles, the inter-particle and 
particle- wall collisions resulting in particle mixing and dispersion, the surface erosion 
due to particle impacts, and particle-turbulence interactions which govern the dis- 
persion and concentration heterogeneities of inertial particles. The dispersed phase 
is usually treated mathematically as a pressureless continuum. Zero-pressure gas 
dynamics was also used for modeling the formation and evolution of traffic jams [3] 
(F. Berthelin, P. Degond, M. Delitala, M. Rascle). 

1.4. Contents of the paper. In this paper we study the problems related with 
the (5-shock in multidimensional system of conservation laws 





pt + V ■ (pU) = 0, 



U t + (U ■ V)U = 0, 



(1.7) 



p t + V ■ ( P F(U)) = 0, (pU) t + V ■ ( P N(U)) = 0, 



(1.8) 
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where F = (Fi, . . . ,F n ) is a given vector field, N = (Ni, . . . , N n ) is a given tensor 
field, N k = (N kl ,...,N kn ), k = l,...,n; Fj, N kj : R» R; p = p(x,t), U = 
[u\{x, t), . . . , u n {x, t)) € R" arc the unknown functions; x = (x\, . . . , x n ) € R n , 
t > 0. System (jl.8p can be rewritten as 

n a n n 

3=1 3 3=1 

The well-known particular cases of this system are zero-pressure gas dynamics in 
the standard form (fl~4|) (here F{U) = U, N(U) = U ® U) and in the relativistic 
form 

Pt + V-{pC{U)) =0, (pU) t +V ■(pU®C(U)) = Q, (1.9) 
(here = C(Z7), iV(J7) = U® C(U)), where Cf?/) = ^=a£ — , Co is the speed 

\/Co + \ U \ 

of light. The relativistic form f) 1 .0[) of zero-pressure gas dynamics was presented 
in [3S]. 

In Sec.[2l we introduce the integral identities (|2.2[) which constitute Definition ^. II 
of i5-shocks for system (jl.8p . Next, using this definition, the Rankine-Hugoniot 
conditions (|2.6p for curvilinear <5-shocks are derived. 

In Sec. [3l geometric and physical aspects of <5-shocks in system (|1.8j) arc studied. 
It is well-known that if U € L°° is a generalized solution of the Cauchy problem com- 
pactly supported with respect to x, then the integral J Rn U(x,t)dx is independent 
of time. For <5-shock wave type solutions this fact does not hold. Nevertheless, by 
Theorems l3.ll "generalized" analogs of these conservation laws are derived. We prove 
that the "mass" and "momentum" transportation processes between the volume out- 
side the moving <5-shock front T t and the front T t are going on. Moreover, we derive 
the (5-shock balance relations (|3.3[) which show that the total "mass" M(t) + m(t) 
and "momentum" P(t)+p(t) are independent of time, where M(t), P(t) are "mass" 
and "momentum" of the domain outside the wave front, and m(t), p(t) are "mass" 
and "momentum" of the wave front T t . 

In SeclU we study the case of zero-pressure gas dynamics. The Rankine-Hugoniot 
conditions (|4.ip for zero-pressure gas dynamics (|1.4p and the Rankine-Hugoniot 
conditions (|4.2p for its relativistic form (|1.9p are particular cases of (|2.6p . For zero- 
pressure gas dynamics system (|1.4p . "mass" and "momentum" have a sense of real 
mass and momentum. In this case the mass transportation process described by 
Theorem 13.11 is the mass concentration process on the moving front Tt (see Theo- 
rem 23]) • According to Theorem 14.21 for zero-pressure gas dynamics energy of the 
volume outside the wave front and the total energy are nonincreasing quantities. We 
consider the possibility of the effects of kinematic self- gravitation and dimensional 
bifurcations of 5 -shock. 

In this section we also consider a spherically symmetric case of zero-pressure gas 
dynamics (|4.3p and present the Rankine-Hugoniot conditions (|4.5p for ^-shocks in 
system (|4.3p . Recall that a spherically symmetric case of the gas dynamics admits a 
solution which describes the heavy shock. This solution related with the investigation 
of atomic bomb explosion was found by L. I. Sedov, J. von Neumann, and G. I. Taylor 
(see [47l 6.16.]). It seems natural that a <5-shock type solution of system (|4.3p can 
model a super explosion with a growing amplitude of the wave front. 

In Appendix [3] some auxiliary facts are given. In particular, we give results 
related with moving surfaces and distributions defined on these surfaces, and prove 
the surface transport theorems. 
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In the author's opinion, the multidimensional system of conservation laws (|1.8|) 
and its generalizations can be used in physical models which can be treated mathe- 
matically as a pressureless continuum, e.g., dusty gas, granular media (for example, 
see [25], [26], [57j, [33]- [35], [2], [32], [20] and Subsec. [L3]). 



2. (5-SHOCK type solutions and the Rankine-Hugoniot conditions 

2.1. (5-Shock type solutions. Let F = {(x,t) : S(x,t) = 0} be a hypersurface 
of codimension 1 in the upper half-space {(x,t) : x G R n , t G [0,oo)} C R n+1 , 
S E C°°(R" x [0, oo)), VS(x, t)\ s=Q for any fixed t, where V = 
Let T t = {a; : S(x,t) = 0} be a moving surface in R". Denote by v the unit space 
normal to the surface T t pointing (in the positive direction) from = {x G R™ : 
S(x,t) < 0} to Q,f = {x G R" : S(x,t) > 0} such that Vj = y^y, j = l,...,n. 
The direction of the vector v coincides with the direction in which the function S 
increases, i.e., inward the domain Qf. Denote by — G = y^y the velocity (along the 



normal v) of the moving wave front L t (see Appendix lA.ljl 
For system (|1.8[) we consider the S-shock type initial data 

(U°(x),p (x);U$(x),xer ), where p°(x) = f{x) + e°(x)S(T ), (2.1) 

U° G L°°(R ra ;R"), p G L°°(R";R), e° G C(T ), T = {x : S°(x) = 0} is the 
initial position of the <5-shock front, VS' (x)| s0 _ ^ 0, U$(x), x G Fo is the initial 
velocity of the 5-shock, S(To) (= S(S )) is the Dirac delta function concentrated on 
the surface To- The facts related to distributions defined on surfaces can be found 
in Appendix IA.21 

Let us introduce a definition of a S-shock wave type solution for system (|1.4I) . 



Definition 2.1. Distributions (U, p) and a hypersurface T, where p(x,t) has the 
form of the sum 

p{x,t) = p(x,t) +e{x,t)6(T), 

and U G L°°(R n x (0,oo);R Tl ), p <E L°° (W x (0, oo); R), e G C(T), is called a S-shock 
wave type solution of the Cauchy problem (|1.8|) , (|2.ip if the integral identities 

<5</? dp(x,t) 



J J p(ft + F(U) ■ Vi^) dxdt + J 



e- 



r ^ V1 + G2 

pP{x)ip(x,0)dx+ I e°{x)ip(x,0)dp,(x) = 0, 
(5</? d/x(a;, i) 



r„ 



p(Ucp t + N(U) ■ Vip^j dx dt+ eU s 



(2.2) 



St VTTcP 



+ J U°(x)pP(x)(p(x,0)dx + j e°(x)U$(x)ip(x,0)dp(x) = 0, 

hold for all <p G T)(R n x [0, oo)), where J f(x)dx denotes the improper integral 
f Mn f(x) dx; 

is the 5-shock velocity, —G = y^y , is the ^-derivative with respect to the time 
variable (|A.5j) . 

It is easy to verify that for n = 1 Definition 12.11 coincides with the definition 
of (5-shocks for one-dimensional zero-pressure gas dynamics (|1.3[) introduced in [TU 
Definition 1.2.1. 
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Let S° be a given smooth function. Denote by 17q = {x : S°(x) < 0} and 
Qq = {x : S°(x) > 0} the domains on the one side and on the other side of the 
hypersurfacc r — {x : S°(x) = 0j. In order to study the (5-shock front-problem, 
i.e., to describe the propagation of a singular front T starting from the initial position 
Tq, we need to solve the Cauchy problem for system (|1.8|) with the initial data 



(p°,U°,U$), where p°(x) = p+°(x) + [p°(x)]H(-T ) + e°(x)6(T ), 

U°[x) = U+°(x) + [U°(x)]H{-T ). [ > 

Here [J7°(x)] = U~°(x) — U +0 (x) is a jump of the function U° across the discontinuity 
hypersurfacc r ; U° = U 0+ , p° = p 0+ if x G n+ , and U° = U°- = U 0+ + [U% 
p° = = p a+ + [p°] if x G £1^; e° and p are given functions, U 0± are given 
vectors; H{— Tq) (= H(—S°))is the Hcaviside function defined on the surface Tq, 
H(-T q ) = 1 if S°(x) < 0, h{-T ) = if S°(x) > 0. Wc assume that for the initial 
data (12. 4p the geometric entropy condition 

U 0+ (x)-v°L <U2(x)-v°L <U°-{x)-u°\ (2.5) 



?"( 

|VS»(x)| 

5 (a;) < 0} to !!+ = {« R fl : S°(aO > 0}. 



holds, where v° = lygo^i is the unit space normal of Tq oriented from Q 



2.2. Rankine Hugoniot conditions. Using Definition 12.11 we derive the 5-shock 
Rankine-Hugoniot conditions for system (|1.8[) . 

Theorem 2.1. Let us assume that Q C R™ x (0, oo) is a region cut by a smooth 
hypersurface T = {(x,t) : S(x,t) = 0} into a left- and right-hand parts flzp. Let 
(U,p), r be a 5-shock wave type solution of system (|1.8| (in the sense of Defini- 
tion al.!^ , and suppose that {U,p) is smooth in fl± and has one-sided limits U , fr 1 
on r. Then the Rankinc-Hugoniot conditions for the 5-shock 

Se 

1 -+Vr t -(eU s ) = ([pF(U)],[p])-n, 

stem ( 2 - 6 ) 

±^l + V rt .(eU s ®U s ) = ([pN(U)],[ P U])-n, 

hold on the discontinuity hypersurface T, where n = (i>, — G) = — fcA— is the space- 
time normal to the surface T, V( Xit ) = (V, Jj), [f(U, p)] — f(U~, p~) — f(U + , p + ) 
is a jump of the function f(U,p) across the discontinuity hypersurface T, 4i is the 
S-derivative (| A.5|) with respect to t, and the tangent gradient Vr t = fjf^) ■ ■ ■ 1 s§~^j 

to the surface T t is defined by (|A.5[) . (|A.6|) . The equivalent forms of (|2.6|) are the 
following: 



or 



Se 

^+Vr t -(eU 5 ) = ([ P F(U)}-[p}U s )-v, 

S -^+\7 rt .(eUs®Us) = {[pN(U)]-[pU}U s )- 

^-2/CGe = ([pF(U)] - [p]U e ) ■ v, 

6(eU s ) 



(2.7) 



(2.8) 



2KGeU & = {[ P N{U)]-\pU]U 5 )-v, 

where IC is the mean curvature (|A.7[) of the moving wave front T t . 

Proof. For any test function ip € V(Q) we have <p(x,t) = for (x,t) ^ G, G C !!. 
Selecting the test function <p(x, t) with compact support in f2±, wc deduce from (|2.2|) 
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that (|1.8p hold in f2±, respectively. Now, if the test function <p(x,t) has the support 
in f2, then 



p hp t + F(U) ■ V<^J <±c 

= / jo(^ t + F(U) ■ V(p) dxdt+ / p(ip t + F(C7) • Vy>) dx dt. 

Using the integrating-by-parts formula, we obtain 

/ p(<p t + F(U) ■ V<p) dxdt = - (p t + V • (pF)) <p(x, t) dx dt 

Jn±nG ^ ' Jn±nG ^ ' 

T I (J^% + ^ U±) 'J S ) <p{x, t) dp(x, t) - [ A*M*, 0) dx, 

JrnG v I V (;c ^)i> | I V (a; ^)i) | t Jn ± nGnw 

where dp(x, t) is the surface measure on T. Next, adding the latter relations and 
taking into account that pt + V ■ (pF) — 0, (a;, t) £ we have 



p(y t +F(U)-Vip} dxdt + J fi{x)ip(x,ti) 



dx 



(-[p]G+[pF(U)]-,)^x,t)^^=. (2.9) 
Now, using the second integrating- by-parts formula in (|A.13[) . one can see that 

where the adjoint operator ^ is defined in (|A.14[) . Thus 



r St VTTG* 



o 



Adding (J2J]) and (|2TT0|) . we derive 

j v (-[p]G+ [pF(U)} ■ v - | - V rt • (cGi/)) <p(x, t) *&?L = 0, 

for all <p(x,t) £ T>(yi). Taking into account formula (|2.3[) for the (5-shock velocity, 
one can see that the last relation implies the first relation in (|2.6p . 

In the same way as above, we obtain the second relation in (|2.6[) . 

In view of (|2.3p and (|A.14[) . the Rankine-Hugoniot conditions (|2.6[) can be rewrit- 
ten as (|2.7[) . Since, according to the proof of Lemma [A. II (see formula (j A. 14|) ^> . 

Vr t • (eU s ) = -2/CGe, Vr t • (eU s <g> E/a) = -2KGeU 6 , (2.11) 

the Rankine-Hugoniot conditions (|2.7p can also be rewritten in the form (|2.8[) , where 
/C is the mean curvature (|A.7[) of the surface T t . □ 

The right-hand sides of the first and second equations in (|2.6[) or (|2.7[) are called 
the Rankine-Hugoniot deficits in p and pC/, respectively. 

Remark 2.1. (a) The Rankine-Hugoniot conditions (|2.6[) constitute a system of 
second-order PDEs. According to this fact, to solve the Cauchy problem for sys- 
tem (|1.8[) . we use the initial data (|2.ip which contain the initial velocity Uj-(x) 
of a (5-shock. It is similar to the fact that in the measure-valued solution ap- 
proach [3D], [3T], [1DJ (see Subscc. |1.2|) the velocity U is determined on the dis- 
continuity surface. 
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It remains to note that according to our Definition ^, ll a <5-shock wave type solution 
is a pair of distributions (U,p) unlike the Definition (|1.5[) . where p(x,t) is a measure 
and U{x,t) is understood as a measurable vector-function which is defined p(x,t) 
a.e.. 

(b) For system (|1.4j) the Rankinc-Hugoniot conditions (|2.6|) have the form (|4.1|) . 
The Rankinc-Hugoniot conditions (|4.1|) arc analogous to the Rankinc-Hugoniot con- 
ditions 



dX 
dt 
dw 
~dt 

d(wU s ) 



U 5 (s,t), 

([pU],[p])-(N x ,N t ), (2.12) 

([pU®U},[pU})-(Nx,N t ), 

in the measure-valued solution approach [3D], [3T], [JDJ (see Subsec. II. 2p . where 
(Nx,N t ) is the space-time normal to the (5-shock front. 

3. Geometrical aspects of 5-shocks: volume balance relations 

It is well known that HUE L°°(R x (0, oo); R m ) is a generalized solution of the 
Cauchy problem (jl.l[) compactly supported with respect to x, then the integral of 
the solution on the whole space 

U{x,t)dx = J U°(x)dx, t>0 (3.1) 

is independent of time. These integrals can express the conservation laws of total 
area, mass, momentum, energy, etc. For a <5-shock wave type solution the classical 
conservation laws (|3.ip do not hold. However, there is a "generalized" analog of 
conservation laws (|3.1[) . In the one-dimensional case these "generalized" analogs 
were derived in [T], [3D], [H]. Now we derive multidimensional generalization of 
these laws. 

Let us assume that a moving surface T t = {x : S(x, t) = 0} permanently separates 
R™ into two parts Slf ={i€l": ±S(x, t) > 0}. Denote fi± = {x £ R™ : ±S°(x) > 
0}. Let (U, p) be compactly supported with respect to x. Denote by 



and 



M(t)= / p(x,t)dx M(0)= / p°(x)dx, 

P{t)= I p(x,t)U(x,t)dx, P(0)= I p°{x)U°{x)dx, 

m(t) = / e(x,t)dp(x), m(0) = / e°(x)dp(x), 

p(t)= [ e(x,t)U s (x,t)dn(x), p{0) = [ e°(x)U$(x)dp(x), 

"masses" and "momentums" of the domains ft^ U ttf , f2^~ U and the "masses" 
and "momentums" of the moving wave front T t , T Q , respectively, where dp,(x) is 
the surface measure on T t . The quantities M(t) and P(t) can be interpreted as the 
volumes under the graphs y = p(x, t) and Y = p(x, t)U(x, t), x £ fl^ U = {x £ 
R" : S(x,t) ^ 0}. 

Theorem 3.1. Let (U,p) and the discontinuity hypesurface T — {(x,t) : S(x,t) = 
0} be a S-shock wave type solution (in the sense of Definition 12. ip of the Cauchy 
problem (|1.8|) . (|2.1|) . compactly supported with respect to x, where p(x,t) = p(x,t) + 
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e(x,t)S(T). Suppose that (U,p) is smooth in il± and has one-sided limits U ± , p^ 
on r. Then the following "mass" and "momentum" balance relations hold: 

M(t) = -m(t), P(t) = -p(t); (3.2) 

M(t)+m(t) =M(0) + m(0), P(t) +p(t) = P(0) +p(0). (3.3) 
Thus the "mass" and "momentum" transportation processes between the volume 
f2f~ U f2j = {x € K™ : 5*(a;, £) ^ 0} and £/ie moving front T t are going on. Moreover, 
the total "mass" M(t) + m(t) and "momentum" P(t)+p(t) are independent of time. 

Proof. Let us assume that the supports of U (x, t) and p(x, t) with respect to x belong 
to a compact K £ R™ bounded by <9iC. Let Kf = n K. By ^ we denote the 
space normal to T t pointing from f2^~ to ilf . Differentiating M(t) and using the 
volume transport Theorem IA.1( we obtain 

M{t) = f %dx+ j Gpdp(x), 

JK~UK+ 01 JdK^UdK+ 

where G = —m^- Since pf + V ■ (p ± F(U ± )) = 0, a; G Iv* and the vectors U ± 

and functions p 1 * 1 are equal to zero on the surface dK t except L t , applying Gauss's 
divergence theorem, we transform the last relation to the form 

M(t) = - [ V • (p-F(U~)) dx - [ V • (p+F(U+)) dx+ [ G[p] dp(x) 

p-F{U-)-vdp{x)+ I p+F{U + )-v dp,(x)+ I G[p]dp,(x) 

([ P F{U)]-v-[p]G)dp{x). (3.4) 

Using the first Rankinc-Hugoniot condition (|2.7p and taking into account that G = 
Us ■ v, relation (|3.4[) can be rewritten as 

M(t) = - jf {[ P F(U)} - [p]Us) -vd^x) = - jf + Vr t • (eU 5 )) d»(x). 

According to the surface transport Thcorcm lA.2[ we have 

, (ft+ v v(*j))#)' 

Thus the first balance relation in (|3.2p is proved. 

Repeating the proof of the first balance relation in (|3.2j) almost word for word, 
we derive the second balance relation in (13. 2|) . 

To complete the proof of the theorem, it remains to integrate (|3.2[) with respect 
to t. □ 

4. Zero-pressure gas dynamics 

4.1. Rankine Hugoniot conditions. According to (|2.6p and (|2.7[) . for zero- pres- 
sure gas dynamics (jl.4p and (jl.9p the Rankinc-Hugoniot conditions have the form 

^T t -(eU s ) = {\pU]-\p]U 5 )-v, 



St 
and 



St 

'- V rt • (eU 5 <g> C/ 5 ) = ([/>[/ ® E7] - [pE7]l7«) ■ i/ 



3 +Vr t -(ei7 5 ) = ([pC(C/)] - \p]U s ) ■ u, 



St 

5 -^^ + V Tt -{eUs®U s ) - ([ P U®C{U)]-[pU]U s )- V . 



(4.1) 



(4.2) 
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respectively. Here according to (|2.11[) . Vr t ■ (eUg) = — 2/CGe, Vr f ■ (eU$ ® Us) = 
-2KGeU 5 . 

In this case the Rankine-Hugoniot deficits in p and pU are the currents of mass 
and momentum, respectively. 

Spherically symmetric case. It is easy to see that the solution of (jl.4| with 
spherical symmetry p = p(r,t), U = u(r,t)-, where r = \x\, x e W l , satisfies the 
following system of equations 

p t + [pu)r H pu = 0, (pu)t + (pu 2 ) r H pu 2 = 0. (4.3) 

r r 

In this case T = {(x,t) €. M™ x [0,oo) : S(r,t) = 0},r t = {ieM": S(r,t) = 0}; 

VS{r,t) = S r f, \VS(r,t)\ = \S r \; v = jfM ; G = -#y; the 5-shock velocity (LO) 

is represented as Us = vG = — rfrf ! x € ^ is easy to sec that if / = f(r,t) 

then formulas (|A.5[) have the form 

Sf df S t df Sf n 

Now formulas (|2.11[) read 

V Tt -{eU s ) = -e^— -, X7 Tt -{eUs®Us) = e ^ n ~p X ' , (a,i)6r. 

Taking into account the above formula, we observe that the Rankine-Hugoniot 
conditions (14.11) take the form 



St St n — 1 r , o r r , St 

et--^-e r -e- — = \pu]-—- + [p]-—- 

O r J r r | o r | |O r | » 

f-e^ -^tf-e^ +e (^l) 2 !^l = [ pu ^] A. + [pu] A 
V SV / * 5V V SV ' r \S r ) r \S r \ | S r | ' 

for (x, t) G T. 

If S(r,t) = —r + <p(i), the Rankine-Hugoniot conditions (|4. 5[) can be rewritten as 

71—1 

e t + 0(t)e r + e4>{t) = -[pu] + \p]<j>(t), 

2-1 ( 46 ) 
(e^)) t + 0(i)(e0(t)) r + e (^(t)) ^— = -[pu 2 ] + [pu]<fi(t). 

4.2. Physical aspects of <5-shocks: mass, and momentum balance relations. 

In this case p > and t/ can be considered as the gas density and gas velocity, 
respectively. Here "masses" M(t), m(t) and "momentums" P(t), m(t) which were 
introduced in Sec. [3] have the sense of real masses and momentums. 

To solve the Cauchy problem, we assume that for its solution the geometric en- 
tropy condition 

U + (x,t) ■ v\ Tt < U s (x,t) ■ v\ Tt < U~(x,t) ■ u\ rt , (4.7) 

holds, where Us is the velocity (|2.3[) of the (5-shock front r t , U is the velocity behind 
the (5-shock wave front and ahead of it, respectively Condition (|4.7|) implies that all 
characteristics on both sides of the initial discontinuity T t must overlap. For t = 
the condition (|4.7p coincides with (|2.5p . 

Theorem 4.1. In the case of zero-pressure gas dynamics (|1.4p the transportation 
process described by Theorem 13.11 is the mass concentration process on the moving 
front T t : 

M(t) = -rh(t), rh(t)>0, P(t) = -p(t), 

M{t) + m(t) = M(0) + m(0), P(t) + p(t) = P(0) +p(0). 
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Proof. It remains to prove the inequality rh(t) > 0. Since the solution (U, p) of the 
Cauchy problem (|1.4p . (|2.1[) satisfies the entropy condition (|4.7p and p > 0, we 
have for the first relation in (14.11) 



^ + Vr t -(ef/ 5 ) = ([pf/]-MC/ (5 )-H rt 

= (p"(C/- -U s )-u + p+{U s -U + )-u)\ Tt >0. 

This inequality and Theorem 13 . 1 1 imply that m(i) = J r (|f + Vr t • (eUg)) > 

and M(t) < 0. In view of these inequalities, in the case of "zero-pressure gas 
dynamics" mass transportation from the volume Vi^ U Clf to the moving wave front 
r t takes place. □ 

4.3. Energy in zero-pressure gas dynamics (jl.4|) . Let us assume that a moving 
surface T t = {x : S(x,t) = 0} permanently separates R™ into two parts flf = {x 6 
K™ : ±S(x,t) > 0}. Let (U, p) be compactly supported with respect to x. Denote 

by 

W{t) = \( p(x,t)\U(x,t)\ 2 dx, w{t) = \f e{x,t)\U s (x,t)\ 2 dp(x), 

energies of the domain U and of the moving wave front T t , respectively (see 
Sec. [3]). The function W(t) +w(t) is the total energy. 

Theorem 4.2. Let (U,p) and the discontinuity hypesurface T = {(x,t) : S(x,t) = 
0} be a 5-shock wave type solution (in the sense of Definition 12. ip of the Cauchy 
problem (|1.4|) . ()2.1|) . compactly supported with respect to x, where p{x,t) = p(x,t) + 
e(x,t)S(T). Suppose that (U,p) is smooth in Q± and has one-sided limits U , pr 1 
on r. Then energies W(t) and W(t) + w(t) are nonincreasing Quantities: 

J t W(t)<0, j t (w(t)+w(t))<0. (4.9) 

Proof. Let us assume that the supports of U (x, t) and p(x, t) with respect to x belong 
to a compact K € R" bounded by dK. Let Kf = fl^ 1 n K. By v we denote the 
space normal to T t pointing from il^ to Qf . Differentiating W(t) and using the 
volume transport Theorem IA.1[ we obtain 

^-p(x,t)\U{x,t)\ 2 dx 
+ at 

Gp{x,t)\U{x,t)\ 2 dp(x)), (4.10) 

dK~UdK? 

where G = 

Since for x G K system fj 1 .4[) has a smooth solution (p > , U ), this solution also 
satisfies non-conservative form (|1.7|) . One can easily verify that (|1.4[1 and (|1.7|l imply 
that 

( P ± |c/±| 2 ) t + v • ( P ± l^l 2 ^) = 0, xe a-±. 

Next, using the last relation, taking into account that the vectors U ± and functions 
p ± are equal to zero on the surface dKf except L t , and applying Gauss's divergence 
theorem to relation (|4.10j) . we transform it to the form 

W(t) = - [ V- (p~ \U-\ 2 U-)dx- [ S7-{p+\U + \ 2 U+)dx+ [ G[p\U\ 2 ] dp(x) 

JkT JK+ JTt 

p- \U-\ 2 U- -vdp{x)+ j p+\U + \ 2 U+ -v dp.(x)+ [ G[p\U\ 2 ]dp(x) 

Ft JTt JTt 
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\U\ 2 U}- \p\U\ 2 ]U s ) -vdp{x). (4.11) 



We also take into account that G = Us ■ v, G = — . 

Since the solution (U,p) of the Cauchy problem (|1.4|) . (|2.1|) satisfies the entropy 
condition (|4. 7[) and p ± > 0, we have 

([p\U\ 2 U] ~ [p\U\ 2 ]U s ) v 

= { P -\U-\ 2 {TJ- -Us)-v + p + \U + \ 2 {Us-U + )-v)\ Tt >0. (4.12) 

Formulas ([4~TT|) . P~T2")) imply that W(t) < 0, i.e., the first inequality in ([43)1 holds. 

Now we will calculate w(t). Taking into account formula (|A.14[) . due to the surface 
transport Theorem IA. 21 we obtain 

w(t) = i jf (j t (e(x,t)\U s (x,t)\ 2 ) + V rt ■ (e(x, t)|CT 4 (x, t)| 2 C/ 5 )) 



-I f-- 

2 / r , VJt 



^ (-(e(a;,t)|tfa(a;,t)| 2 ) - i)|?7 5 (x,t)| 2 ) 

1 /" ( ( _ S(eu S k) , _fot,5 
2 



/ (E(^^^+^e^)-2/CGe( ; r 1 i)|%(x,t)| 2 )d / i(a;). (4.13) 
•^ r * fe=i 



77«5fc + - 2JCGeu S k = [pu k U ■ v] - [pu k ]U s ■ v, 

St g Jt (4.14) 

— u S k - 2KGeusk = [pU ■ v]u Sk - [p]U s ■ vu Sk , 

where us k (x, t) is the fc-th component of the vector Us, k = 1, . . . , n. Now, subtract- 
ing one equation from the other in (|4.14p . we obtain 
Su 

e-rr- = [p^U ■ v\ - [pu k \Us -v -[pU ■ v]usk + [p]Us • vu &k - (4.15) 
Substituting equations (|4. 15[) into (|4.13[) . one can easily calculate 
1 f n 

w(t) = - / {[pu k U ■ v] - [pu k ]U ■ v)u Sk 

k=l 



-[pU ■ y]\U s {x, t)\ 2 + [p]\U s {x, t)\ 2 U s ■ v) dp{x). 
Taking into account that Us = Gv, we rewrite the above relation in the form 



w(t) = lf (2[p(U ■ v) 2 ]G - 3[ P U ■ v\G 2 + [p]G 3 ) dp(x), (4.16) 



where G — — j^gy ■ 



Adding (|4.11j) and (|4.16|) . we obtain 

Wit) + w(t) = -\l {\P\U\ 2 U ■ v] - [ P \U\ 2 ]Us ■ v 



- 2[p(U ■ vf]G + 3[pU ■ v\G 2 - [p}G 3 ) dp(x). (4.17) 
Let us represent the velocity on the wave front U\r t as the sum of the normal 
component U -v and the component Ut an tangential to the surface T t . Since |t/| 2 |r t = 
[U ■ v) 2 + U 2 an , and G = Us ■ v, one can represent the integrand in (|4.17p as 

([p\U\ 2 U ■ v\ - [ P \U\ 2 ]Us ■ v - 2[p(U ■ v) 2 ]G + 3[pU ■ v\G 2 - [p]G 3 ) 

= p-{Ut an ) 2 {U- . v -U s -v) + p+{U+ n ) 2 {Us -V-U+-V) 

+ p-(U~ -v-Us-vf + p + {Us-v-U- -vf. (4.18) 
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Since the solution (U, p) of the Cauchy problem Q1.4p . (12. ip satisfies the entropy 
condition (|4.7p and p > 0, we deduce that the last expression is non-negative. 
Formulas (|4.17p . (|4. 18|) imply that W(t) + w(t) < 0, i.e., the second inequality in 
(PI) holds. □ 



Remark 4.1. In [10j . (5-shock type solutions of the Riemann problem for one- di- 
mensional zero-pressure gas dynamics were studied as the vanishing pressure limit 
of solutions to the Euler equations for nonisentropic fluids. As stated in [TQ1 p. 142], 
the limit system formally is becomes the system of transportation equations (|1.3p 
with the additional conservation law 

(pE) t + {puE) x = 0. (4.19) 

Next, it is stated p~0l p. 143] that the additional conservation law (|4.19p actually 
yields the entropy inequality 

(pu 2 ) t + {pu% < 0. (4.20) 

in the sense of distributions for the Riemann solutions to Q1.3P . 

According to Theorem 14.21 (for n = 1), the inequality (|4.20p is not the entropy 
inequality, it reflects the fact of energy nonincreasing. 

4.4. Two possible effects. By analyzing the o~-shock Rankine-Hugoniot conditions 
(|4.ip for zero-pressure gas dynamics (|1.4p . one can deduce the possibility of the 
following interesting effects. 

The effect of kinematic self-gravitation. According to (|4.8p . in zero-pressure 
gas dynamics (|1.4p the mass concentration process on the moving discontinuity sur- 
face T t is going on. Moreover, the second (5-shock Rankine-Hugoniot condition in 
(|4.ip is the momentum conservation law. Taking this fact into account and us- 
ing Newton's second law of motion, one can introduce an "effective" gravitational 
potential in a neighborhood of the discontinuity surface and describe the concentra- 
tion process in terms of gravitational interaction. Since in system (|1.4p there is no 
term related with gravitational interaction, this "gravitational effect" is of a purely 
kinematic nature. 

Dimensional bifurcations of o~-shock. It follows from Theorem 14.11 that in 
the n-dimensional zero-pressure gas dynamics (jl.4p the mass transportation process 
from the volume fi^ U f2^~ = {x € M™ : S(x,t) ^ 0} onto the n — 1-dimensional 
moving (S-shock front T t is going on. Let us suppose that in a finite time period t 
the whole initial mass M(0) may be concentrated on T t . Then, according to (|4.ip . 
for t > t, instead of the whole "initial" n-dimensional system of zero-pressure gas 
dynamics (|1.4p we obtain a "surface" (n — 1)- dimensional version of this system 

^ + Vr t • (eU 5 ) = 0, + Vr ( ■ (eU s ® U&) = 0, (4.21) 

dt dt 

where instead of the gas velocity U we have the velocity Us of the moving (5-shock 

front r t , and instead of the gas volume density p we have the surface density of 

the front mass e. Moreover, the quantities Us, e are defined only on the moving 

front r t . Since system (|4.2ip is an n — 1-dimensional analog of system (|1.4p on 

the (n — 1)- dimensional surface T t as on a Riemannian manifold, therefore its 

solution can develop singularities within a finite time period, and the whole mass is 

concentrated on the manifold of dimension n — 2, and so on. Thus, it may happen 

that after a finite number of bifurcations the whole initial mass will be concentrated 

at a singular point. 

A description of the above effect prompts to generalize Definition 12.11 and intro- 
duce a new concept of a multidimensional 5-shock type solution to system (|1.8p as 
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a pair (U, p), where p(x, t) has the form of the sum 

n 

p(x, t) = p(x, t)+J2 ej(x, t)8(rW), (4.22) 
3=1 

U G L°°(R n x (0,oo);R"), p G L°°(R" x (0,oo);R), e G C(rW), T« is a hy- 
persurface of codimension j, S(T^) the Dirac delta function concentrated on the 
hypersurface , j = 1, 2, . . . , n. For this purpose we need to derive special integral 
identities analogous to (|2.2| and develop the theory of such type of solutions. In the 
framework of such type definition one can solve the above problem of dimensional 
bifurcations of <S-shock. 



Appendix A. Some auxiliary facts 

A.l. Moving surfaces of discontinuity. Let us present some results concerning 
moving surfaces from 5.2.], 0], [5]. Let T t be a smooth moving surface of 
codimension 1 in the space R n . Such a surface can be represented locally either 
in the form T t = |.t G R n : S(x,t) = 0}, or in terms of the curvilinear Gaussian 
coordinates s = (s±, . . . , s n -i) on the surface: 

Xj = xj(si, . . . , s n _i, t), s £ l n_1 . 

We also consider the surface T = {(x,t) G M n+1 : S(x,t) = 0} as a submanifold of 
the space-time R™ x R. We shall assume that VS(x, t) | r ^ for all fixed values 
of t, where V = (g§^-, ■ ■ ■ , gf - )- Let v be the unit space normal to the surface T t 
pointing in the positive direction such that J^- = |VS|i/j, j = 1, . . . ,n. 

Let f(x,t) be a function defined on the surface T t for some time interval, and 
denote by 4f the derivative with respect to time as it would be computed by an 
observer moving with the surface. This derivative has the following geometrical 
interpretation. Let M be a point on the surface at the time t = t . Construct the 
normal line to the surface at Mq. At the time t = to + At, At is an infinitesimal, this 
normal meets the surface T t +At at the point M = M(t + At). Then the (5-derivative 
is defined as 

Sf(M ,t ) = Hm fjM) - /(Mp) 

5t At^o At [ ' ' 

If As is the distance between Mo and M, then 

G= Um ^ (A.2) 

is the normal velocity of the moving surface T t and 

^ = Um ^ = lim = Gi>j, j = l,...,n. (A3) 

<5f At^o At At^o At As v ; 

Since it is essential that the ^-derivative is computed on a surface, and S remains 
constant on this surface then ff = 0. Thus we have 

SS dS \ -< SS Sxj dS 9 

3 = 1 J 3=1 

i.e., 

S t = -G\VS\. (A.4) 

From this formula we can see that — G = jy^y can be interpreted as the time 
component of the normal vector. 
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The space-time unit normal to the surface T is given by n = -^==^, where 

v (M) = (v,f). 

If f(x, t) is a function defined only on T, its first order ^-derivatives with respect to 
the time and space variables are defined by the following formulas [23l 5. 2. (15), (16)]: 

Sf def df df Sf def df 9f 

-77 = 77T+(J77-, 7— = a V i~*~i J = 1 ,---> n i ( A - 5 ) 

ot at ov oxj oxj ov 

where / is a smooth extension of / to a neighborhood of T in R™ xt, j = 1, . . . , n, 
and ^ = v ■ V/ is a normal derivative. Thus the gradient tangent to the surface T t 
is defined as 

/ S 6 



Vr <= v - v <-fe'--< *rj> (A - 6) 

where = v{y ■ V) is the gradient along the normal direction to the surface T t . 

Note that the (^-derivatives (|A.5[) depend only on the values of / on T, i.e., if 
/ = on T then j£- and on T, j = l,...,n. Indeed, let {x ,t ) G V. If 

^(x,t)f(xo,k) = then V rt /(a:o,<o) = and &(x ,t ) = 0, where V( X)t ) = (V, ^). 
If V/(xo, to) 7^ then in a neighborhood of the point (xo, to) the surface T t has the 

unit space normal v = an< ^ ^ = — TWT Consequently, Vr t /(xo,£o) = and 
^(x , to) = 0. In the sequel we shall drop tilde from /. 

For a vector A(x,t) = (Ai(x,t), . . . ,A n (x,t)) defined only on T t , we introduce 
the surface (tangent) divergence by the following formula 

n 

3=1 3 

The mean curvature of the surface Tt is defined as 

^ def 1„ 1 A fo, 1„ , . 

A. 2. Distributions defined on a surface. Consider some facts concerning distri- 
butions defined on a surface [531 5.2.], J5TJ ch.III,§l.], 0], [5]. The Heaviside function 
H(S) is introduced by the following definition: 



(H(S), ip(x,t)) = (p{x,t)dxdt, y ip e V(R n x R). 

JS>0 

According to [231 5. 3.(1), (2)], we now introduce the delta function 5(S) on the surface 
F, whose action on a test function ip(x, t) € 2?(R" x R) is given by 

(5(S), <p(x,t)) = ^ ( V {x,t)dn{x)dt= j <p(x,t) %t4= , (A.8) 

where d/i is the surface measure on the corresponding surface. According to |231 
5. 5. Theorem 1.], we have 

— = VJ S{S), = -GS(S). 

Now we introduce the derivative of the delta function d v 8(S) along the space 
normal v by the formula [23j 5.3.(7)] 

(dJ(S), <p)=-(5(S), ^f)=- j f ^dfx{x)dt, V^P(M"xM), (A.9) 
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where ^£ = v . i s the normal derivative of ip. If /(x, t) is a continuous function 
defined on T which is a a restriction of some continuous function defined in a neigh- 
borhood of r in R n x R, then the distribution d v (fS(S)) (the so-called double layer) 
is a functional acting by the rule 

(d v (fS(S)), p) = -(S(S), f^), e P(R n x R). 

According to [23l 5.3.(6)], we have 

™ 8 

5'{S) = J2 Vi^-5(ff) = 2KS(S) + d v 8{S) 



dxi 



and 



— S(S) = -G(2K5(S) + d u S(S)), —5(S) = p j (2K5(S) + d u 5(S)), (A.10) 

where /C is the mean curvature (|A.7[) of the surface T t . 

If f(x,t) is a diffcrcntiablc function, using (|A.5[) . (|A.10|) . one can prove the fol- 
lowing relations [H 12. 6. (15), (16)] 

!(/<*(£)) = (% + g¥-~ 2!CGf)5(S) - Gfd u 8{S). (A.12) 



<9£ ' V 9i <9f 

A. 3. One integrating-by-parts formula. 

Lemma A.l. Suppose that e(x,t) is a compactly supported smooth function defined 
only on the surface T = {(x,i) : S(x,t) = 0}, and e(x,t) is the restriction of some 
smooth function defined in a neighborhood ofT in R" x R, and Tq = {x : S(x,0) = 
0} . Then the following formula for integration by parts holds: 

8<p dfj,(x,t) f S*e dfj,(x,t) f 1Q , 

for any p € 2?(R™ x [0,oo)), where 4r is the adjoint operator defined as 

/C is i/ie mean curvature (|A.7|) o/ i/ie surface T t . 

Proof. With the help of formulas (|A~8|) , fO|) . ([ATTO]) . (fATTT]) . (TA~T2l) . we derive by 
simple calculations 



-(^(e*(S)F(t)), p)-(a,(Ge*(S))F(i), p 



r <Jt a/TTG^ 



(^(S) -eG(2K6(S) + d u S(S)), p) - (e8{S)8{t), <p) 

S(S)ir S -^v k + eGd„S(S), cp 
^ ox k 



fe=i 



/8e 
\8t 



2KGe)6(S), <p\ - (e(x,0)8(S{x,0)), <p(x,0)\, 
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where H{t) is the Heaviside function. Here we use the obvious relation 

fe=l K 

Using the last relation and formula (|A.T|) . we calculate 



3 = 1 J 3=1 



□ 



A. 4. Transport theorems. Here we give the following transport theorems. 



Theorem A.l. ([13 12.8.(3)], g], [6], 0) iei /(a:,*) be a sufficiently smooth func- 
tion defined in a moving solid fl t , and let a moving hypersurface dflt be its boundary. 
Let v be the outward unit space normal to the surface dfl t and W(x, t) be the velocity 
of the point x in Q t . Then the volume transport theorem holds: 

4 / f(x,t)dx = f %dx+ f fW-vdn(x) 
at J fi t Jn t at Jdvtt 

Theorem A. 2. ( [UJ 12-8.(9)]) If e(x,t) is a smooth function defined only on the 
moving surface T t — {x : S(x,t) = 0} (which is the restriction of some smooth 
function defined in a neighborhood ofT t ), then the surface transport theorem holds: 



37 r 6 



dt 



= J r (^~ 2/CGe) dfi(x) = jf (| + Vr t • (eU s )) d^(x), (A.16) 
where Us is the velocity ofT t . 
Proof. Since according to definition (IA.8[) , 

m(t) = / e(x,t) dfj,(x) = (e(x,t)S(S), l) , 

using (|A. 12^ . we obtain 

rh(t) = ^(e{x,t)5(S)), l)^ = ((^ - 2KGe}6(S) - Ged v 6(S), l) ^ 

= {(^-2KGe)s{B), l) x = l t (g-aCGe)^(x). 
To complete the proof of the theorem, it remains to use formulas (|A.14[) and 
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